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ABSTRACT 

This article reviews some effects of disorder in percolation systems even away 
from the critical density p c . For densities below p c , the statistics of large clusters 
defines the animals problem. Its relation to the directed animals problem and the 
Lee- Yang edge singularity problem is described. Rare compact clusters give rise 
to Griffiths singuraties in the free energy of diluted ferromagnets, and lead to a 
very slow relaxation of magnetization. In biassed diffusion on percolation clusters, 
trapping in dead-end branches leads to asymptotic drift velocity becoming zero 
for strong bias, and very slow relaxation of velocity near the critical bias field. 

1 Introduction 

Let me start by thanking Professor H. R. Krishnamurthy and other members of 
the organizing committee for inviting me to this meeting to felicitate Professor 
Narendra Kumar on his 60th birthday, and giving me an opportunity to pay my 
tribute to him. 

Over years, I have always enjoyed discussing various questions with Narendra. 
The large spectrum of his interests, and his spirit of enquiry, and insights to 
disentangle the essential problem from confusing camouflage have been a source 
of admiration for me. Thus I am really very happy to come here, and express my 
respect for him, and join the other speakers in wishing him many more years of 
happy questing. 

While I have shared with Kumar a common interest in understanding disor- 
dered systems, my own work has been largely in classical statistical mechanics 
(H — 0), while Kumar's contributions to quantum problems have been discussed 
by several speakers here. Even so, I can legitimately claim to be one of Kumar's 
coworkers. Our only paper together was on the behavior of ±J spin-glass on 
a Bethe lattice, and was presented at the DAE Solid State Physics Symposium 
meeting at Madras ( now Chennai) in December, 1979. It turned out that neither 
Narendra nor I could attend the meeting, so the paper was presented by our host, 



Profesor Krishnamurthy. The abstract of the paper appeared in the proceedings 
of the symposium. Unfortunately, both of us got enmeshed in other problems, 
and did not ever write up this work for publication in regular journals. 

Professor Krishnamurthy has suggested that it would be better if I discuss 
generally the different effects of disorder in classical statistical mechanics, and 
not narrowly focus on the topic of my latest preprint. I was rather relieved by 
this, as my paper-writing has been even slower than my usual (slow) rate in the 
recent past. So, in this talk, I will not describe any new results, but discuss 
some examples of interesting effects caused by disorder in a classical percolation 
system. The selection of topics was determined primarily by my familiarity rather 
than any other reasons. I will only outline the main results. The interested reader 
will have to go back to the cited literature for details. 

Consider a random mixture of equal-sized conducting and insulating beads in 
a box. It is easy to see that as the fraction of conducting balls is varied from zero 
to one, the bulk mixture undergoes a transition from insulating to conducting. 
This is the basic percolation transition, and has been studied a lot as a simple 
geometrical model of phase transitions. As in the case of thermal critical phase 
transitions, various physical quantities have singular behavior near the phase 
transition, characterized by "critical exponents". Scaling theory, simulations and 
theoretical techniques such as renormalization group etc. have been used very 
successfully to understand the behavior of percolation systems at its critical point. 
In particular, in two dimensions, all the critical exponents of percolation are 
known from conformal field theory. Several good reviews of this subject are now 
available [|IJ. 

There is more to percolation theory than the critical exponents. Of course, 
an experimental disordered system may often be modelled by site- or bond- per- 
colation. In general, one is much more likely to find it not near its percolation 
threshold. Critical exponents of the percolation theory are not of much use in de- 
scribing these systems. It turns out that percolation systems, as classical models 
of disordered media show many interesting features, even away from the critical 
point. I will discuss some examples of these in the following. 

2 Off-critical Exponents in Percolation 

Let P n {p) be be probability that the cluster of connected sites containing the 
origin has exactly n sites. Then, for all p < p c , for large n, 

P n (p) ~ An~ 9 exp[-B(p)n] (1) 

where B(p) is a p-dependent function that goes to zero as p tends to p c . The 
exponent 9 is independent of p, and depends only on the dimension of space. This 
law is valid for n much greater than the typical cluster size n*(p). For n <C n*(p), 
one gets a different exponent P n {p) ~ n~ T , with r ^ 9. As p tends to p c , n*(p) 
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diverges. Here r is a critical exponent of percolation theory, but 9 is an off-critical 
exponent. 

In the limit p — > 0, all clusters of n sites have the equal weight ( p n ). Let 
A n be the number of different clusters of n sites possible that contain the origin. 
(These are called rooted animals: same cluster with different positions of origin 
are counted as distinct @.) Then one expects that for large n, 

A n ~ KX n n- e (2) 

The exponent 9 in this equation is same as that defined by Eq.(|l|), because for 
small p, P n {p) ~ A n p n . One can also define the average linear size of an animal 
of N sites. This grows as n u , where the exponent v is related to the exponent 9 
defined above by the relation 

9=(d- 2)v (3) 

where d is the dimension of space. The above equation is valid for 1 < d < 
8. For d > 8, the exponents 9 and v stick to their mean-field values 3/2 and 
1/4 respectively. Eq.(3) has the form of a hyperscaling relation except that 
[d — 2) appears here instead of d. This is understood as being due to a hidden 
supersymmetry in the problem y], which makes the problem of determining the 
number of animals in ^-dimensions related to the problem of Lee- Yang edge 
singularity in (d — 2) dimension. 

The Lee- Yang description of the mathematical mechanism of phase transitions 
is well-known j|] . For a hard-core lattice gas in a finite volume, with possible addi- 
tional attractive short-ranged interactions, the grand-canonical partition function 
is a finite degree polynomial with positive coefficients in the chemical activity z. 
The zeroes of this polynomial often ( not always) lie on lines in the complex- z 
plane. As the temperature is varied, the coefficients of the polynomial change, 
and the zeroes move. If at some temperature, the zeroes come arbitrarily close 
to the real axis as the size of the system is increased, the free energy per site 
becomes a non-analytic function of z, signalling the onset of a phase-transition. 

The density of zeroes along such a line of zeroes near its end point shows a 
power-law dependence on the distance from the endpoint. We define the Lee- 
Yang edge singularity exponent a by the relation that the density varies as e u at 
distance e from the endpoint. It turns out that a is independent of temperature 
for all temperatures T above the criticial temperature T c , and depends only on 
the dimension of the system. It was shown by Parisi and Sourlas |3j that the 
animal exponent 9(d) in ^-dimensions is related to the Lee- Yang edge singularity 
exponent in d — 2 dimensions 

9(d) = a(d - 2) + 1 (4) 

If we allow only neighbors in the "forward direction" , ( say along the direc- 
tion of increasing coordinates on a hypercubical lattice), we get animals with a 



3 



directional constraint. It turns out [0] that critical exponents of directed animals 
in d dimensions turn out to be related to those of undirected animals in (d + 1) 
dimensions. 

e dtr (d) = e(d + i) (5) 

and the transverse size exponent v^^ ir for directed animals in d dimensions is 
the same as the (only one) size exponent v for undirected animals in [d + 1) 
dimensions 

v±,dir{d) = v{d+ 1). (6) 

The exponent a is easily shown to take the values a = — 1 for d = ( a 
point), and o = —1/2 for d = 1. For d = 2, one can use the exact solution 
of hard hexagon lattice gas by Baxter to show that a = —1/6 for (i = 2 |. 
This then shows that the exact values of the exponent 9 for undirected animal in 
dimensions 1, 2, 3, 4 are —1, 0, ~, | respectively. The corresponding values of the 
size exponent v are 1, ~, ^ in dimensions 1,3,4 respectively. In two dimensions, 
the exponent v is not determined by scaling relations given above. The upper 
critical dimension for the animals problem is 8, and for all d > 8, we get 9 — |, 
and v — j. The exponents for directed animals are easily determined from the 
scaling relations given above. 

We have seen that there is a fairly good understanding of the off-critical 
"below p c " exponents of percolation. One can also define off-critical exponents 
in the super-critical regime of percolation theory. It was shown by Kunz and 
Suillard that for all p > p c , the probability that the origin belongs to a finite 
cluster of n sites varies as exp(—b(p)n~d~) , for sufficiently large s. This is easy to 
understand: to get a finite cluster of n sites, we need to disconnect a it from the 
infinite cluster. This needs perimeter bonds order ( n~) bonds, which gives the 
result. More accurately, there is apower-law prefactor multiplying the exponential 
term, and the probability for a finite cluster of n sites varies as 

P n (p) ~ Kn~ 6 exp(—b(p)n~) (7) 

where b(p) is a function of p, and K is a constant. The exponent 8' can be 
calculated exactly || using the fact that such clusters are roughly compact, with 
linear size varying as n^, and the fluctuations in the (d — l)-dimensional roughly 
spherical surface can be described in terms of normal modes of vibration of the 
surface. Lubensky and McKane using field theory techniques showed that the 
exponent 9' takes the values |, — ~, |, — — 1| for d = 2, 3, 4, 5, 6 respectively. 
The non-monotonic behavior of 9' as a function of d comes from the fact that the 
integrals over normal modes in odd and even dimensions coming in the theory 
have different behaviors. 

Note that this is one of the few cases where non-trivial values of exponents 
can be exactly calculated in many dimensions greater than 2. 
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3 Relaxation in Disordered Ferromagnets 



Consider now the case when the atom at each occupied site in a percolation net- 
work carries a magnetic moment, and there is a nearest neighbor ferromagetic 
interaction J between the magnetic atoms. The "unoccupied" sites may be va- 
cancies, or occupied by nonmagnetic atoms. Let p denote the concentration of 
magnetic atoms. Then if p < p c , the percolation threshold, then the system 
breaks up into mutually-disconnected clusters of magnetically coupled spins. In 
such a system no long-range spontaneous magnetization is possible at any non- 
zero temperature. For p > p c , there is an infinite connected cluster, and at 
sufficiently low temperatures, spontaneous magnetization exists. The transition 
temperature T curie (p) depends on p, and goes to zero, as p is decreased from 1 to 
p c . The ferromagnetic phase is denoted by F in Fig. 1. 

For p < p c , there is no spontaneous magnetization. Let f(T, h) be the 
disorder-averaged free energy per site of this system at a temperature T in a 
magnetic field h. However, it was shown by Griffiths that f(T,h) is a non- 
analytic function of the magnetic field at h = for all T < T curie (p =1). While 
the partial derivatives ^ r /(T, h) exist for all positive integers n, and are finite, 
the Taylor series for f(T, h) in powers of h does not converge for any T below the 
T cur i e {p = 1), the Curie-temperature of the "pure" system. Thus, in the entire 
region marked G ( for the Griffiths phase) in fig. 1, there is no spontaneous 
magnetization, but the free energy per site f(T, h) is a non- analytic function of 



While the nonanalyticity in f(T, h) as a function of h is of the nature of 
an essential singularity, and is difficult to verify in experimental systems, the 
rare large clusters responsible for it have a much more pronounced effect on the 
rate of relaxation to equilibrium in disordered systems. Consider the decay of 
magnetization in such a system. We assume that the system is coupled to heat 
bath at temperatute T < T curie (p = 1), and relaxes by single-spin -flip Glauber 
dynamics. At time t = 0, the system is prepared in a state with all spins up. We 
monitor the bulk magnetization M(t) at large times t. Since the spins in different 
clusters do not interact with each other, we can write M(t) as a weighted sum 
over different cluster configurations C. 



For a given finite cluster C of n-sites, one can determine the average magnetization 
at a site in the cluster by explicit diagonalization of the of matrix of size 2" x 2 n . 
This is non-trivial, except for very small n. Fortunately, the behavior of M(t) for 
large t can be determined by simple qualitative arguments. 

In the summation Eq.(§), the leading behavior of each term is exponential in 
time, with the decay rate 7c depends on the cluster C. Thus, we may write 



h. 



M(t) = ^Prob(C) < S(t) > c 



(8) 



c 



< S(t) > c 



exp(-j c t) 



(9) 
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and hence 

M(t) ~ Prob(C)exj9(-7c t) (10) 

c 

At large times, clusters with smallest decay rates contribute most. The slowest 
relaxing clusters are those where all sites within a disc of radius R are occupied. 
The density of such clusters varies as exp(—AR d ), where the constant A depends 
on p. If the magnetization of such a cluster has to flip, it would need to create a 
domain wall of energy ~ cxi? d_1 . The rate of such activated process will decrease 
for large R as exp( — aR^ 1 /T). It is then straight forward to put these estimates 
in the Eq.(g), and deduce that 

M(t) ~ exp(-K(logt)&), for large t. (11) 

For p > p c , in the ferromagnetic phase (F in fig. 1) the probability of a large 
finite cluster of radius R, varies as exp(— ai? d_1 ), and not as exp(—R d ). Then the 
steepest descent calculation shows that in the ferromagnetic phase F, the mag- 
netization at long times decays as a power law M(t) ~ t~ c , where the exponent 
c depends on both p and T. 

The argument outlined above was first presented in pi. The argument has 



been refined ||TT| , [TI|. Unfortunately, neither actual experimental data on disor- 
dered ferromagnets, nor results of numerical simulations |]TB| show a clear evidence 
of such a exp[—(logt) x ] behavior. Presumably the time scale beyond which the 
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contribution of rare clusters will dominate is larger than experimentally accessible 
time-scales. The experimental data seems to fit better a stretched exponential 
exp(—t x ). It seems that a more careful argument, that gives not only the correct 
asymptotic behavior at longest times, but also at intermediate times is needed. 



4 Biassed Diffusion on Percolation Networks 

Consider the motion of a single diffusing particle on a percolation network, say 
in two dimensions, with density of occupied sites being p. We assume that the 
particle can move only on the occupied sites of the lattice. Then if p is less than 
the critical probability p c , the particle is localized. For p > p c , if the particle 
starts on the infinite cluster, its mean square deviation from the initial position 
grows linearly with time, < R 2 >~ D(p)t, where the diffusion constant D(p) 
depends on p and tends to zero as p tends to p c . This problem of (unbiassed) 
diffusion on percolation clusters has been studied a lot fl4j . 



If there is a larger probability of displacement in some direction, due to an 
imposed field, we have biassed diffusion. We shall model it by assuming that at 
any time the diffusing particle attempts to take a step in the up, right, down, left 
directions with probabilities (1 — B)/A, 1/4, (1 + B)/4 and 1/4 respectively. The 
step is actually taken if the intended destination site is occupied. If the biassing 
field B is small, we have a non-zero value of average displacement per step, and 
this gives rise to mean displacement in time t in the direction of the field growing 
linearly with t, and the mean velocity in time time t tends to a constant 

— * 

{Joo = Lim^oo < R t > jt (12) 

This asymptotic velocity is proportional to B for small B. 

If p is near 1, most of the sites are occupied, and at large length scales, the 
medium looks homogenous. One then expects that so long as p > p c , we expect 
the same behavior as in the system without disorder ( p = 1). 

Bottger and Bryksin realized that this is not so [|IJ|]. They argued that the 
mean velocity must tend to zero as B tends to 1, because of the possibility of 
trapping in dead-ends. We argued in [|T(J, that for any p < 1, there exists a 
critical bias B c such that the asymptotic drift velocity v is exactly zero for all 
B > B c . This is easily seen: during its motion, the particle may get trapped in 
dead-end branches for long times, as it has to move against the field to get out of 
the trap. For a trap of depth £ the potential barrier to cross increases with £, and 
the trapping time varies as (jzjjY- The density p(£) of traps of depth £ varies as 
exp(—£/C,). Hence the average trapping time per step along the backbone is 

OO -I I p 

YsPWij^y (is) 
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Figure 2: 



where £ is the p-dependent percolation correlation length of the system. This 
summation converges only for B < B c = tanh(^). For B = B c — e, with e > 0, 
this summation varies as 1/e, and the mean velocity, which varies inversely as 
the the mean trapping time varies proportional to e. 

For B > B c , the asymptotic velocity is zero, but the mean displacement 
of the particle < R t > increases as t a with a < 1. The exponent a depends 
continuously on B, and is easily obtained by using the condition that the density 
of traps with trapping time geater than or equal to t varies as t~ a , where 

a=(|)M^|) (14) 

In a time t, the particle can, on the average, only travel a distance t a before 
it encounters a trap with trapping time bigger than t, and gets stuck there. 
Eventually, it will exit from this trap, only to get stuck in other traps, some 
with an even larger trapping times. Thus, if we examine all particles at some 
large time t, typically they would be stuck in, or just emerging from a trap with 
trapping time of order t. 

For B = B c , the average trapping time per step encountered by the walker is 
given by (|I3D, except that summation over £ is cutoff at a value i ma x, where £ max 
is the typical value of the trapping time of the deepest trap encountered by a 
walker upto time t. It is easy to see that £ ma x varies as logt for large t, and hence 
the average velocity of the particle up to time t varies as y^-r. Thus the average 
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velocity decreases very slowly to zero at B = B c . This very slow relaxation has 
been checked in large-scale simulations of this problem |17 . 



If B = B c — e, then this slow decay of the velocity continues from its initial 
value of order 1, to the final value which is of order e. As the initial decay of 
velocity would be nearly same as that for B = B c , we see that typical relaxation 
time t(B) for the average velocity in an ensemble of non-interacting particles to 
reach the steady state value varies as 

t(B) ~ exp{^-^) (15) 

where A is a constant. Note that this relaxation is the overall relaxation time 
for a macroscopic observable ( average current- density for an ensemble of non- 
interacting particle). We are not discussing the largest relaxation time, as that 
is infinite for all B, however small. 

Thus, biassed diffusion of noninteracting particles on a percolation network 
provides a very simple model where a fast rise of relaxation time near a dynamical 
phase-transition ( " the Vogel-Fulcher law" of glassy dynamics) can be seen. 

We see that disorder effects both the static and dynamical properties of the 
system in a very significant way The effect on non-equilibrium properties like 
response functions is much more pronounced. I hope that further work will lead 
to a better understanding of these systems. 
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Captions to figures 

Fig. 1: The phase digram of a diluted magnet. The ferromagentic, paramagnetic 
and Griffiths phases are denoted by PM, F, and G respectively. 

Fig. 2: Schematic representation of the percolation cluster with density above the 
critical threshold. The heavy lines denote the dead-end branches. 

Fig. 3: The average velocity v as a function of the biassing field B. 
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